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I. INTRODUCTION 



Electron scattering from nuclei is a common and powerful tool for studying the structure 
of nuclei. The method relies on our superior understanding of Quantum Electrodynamics 
(QED) and the relative ease by which QED may be applied to a variety of processes, at 
least in the one-photon-exchange approximation. In inclusive (e, e') electron scattering all 
nuclear-structure information is contained in two dynamical quantities: the longitudinal and 
transverse response functions. The longitudinal response is sensitive to the distribution of 
charge in the nucleus while the transverse response samples the distribution of currents and 
magnetization. Measurement of these quantities in the quasielastic region is expected to be 
particularly clean as the reactive content of the reaction is dominated by quasifree proton 
knockout. If so, "reduced" longitudinal and transverse response functions, obtained from 
the full nuclear responses by dividing out the corresponding single-nucleon form factor, 
should be equal. Yet a quenching of the longitudinal response relative to the transverse 
one of 14% in 4 He and 50% in 208 Pb has been reported from a quasielastic (e, e') electron- 
scattering measurement A similar, in fact even larger 20-40%, quenching in 4 He has also 
been reported in the semi-exclusive (e, e'p) reaction at quasielastic kinematics 0. In order 
to explain the longitudinal/transverse (L/T) discrepancy a variety of explanations have 
been proposed. These include medium-modifications to vacuum polarization ||, nucleon 
"swelling [[|, and Brown- Rho scaling ||. It is fair to say, however, that the L/T anomaly 
remains an unsolved problem in nuclear physics. 

The appeal of the (e, e'p) reaction is due to the perceived sensitivity of the process to 
the nucleon momentum distribution. Interest in this reaction has stimulated a tremendous 
amount of experimental work at electron facilities such as NIKHEF, MIT/Bates, and Saclay, 
who have championed this effort for several decades. While it is undeniable that this reaction 
involves the best understood theory in all of physics (QED) many uncertainties remain due to 
the strongly-interacting character of the many-body system. It is hoped that with the advent 
of modern electron-scattering facilities, such as the Thomas Jefferson National Accelerator 
Facility (JLab) and Mainz, some of the remaining open questions will be answered. Indeed, 
in an attempt to elucidate the physics of the L/T anomaly discussed earlier, a systematic 
study of the longitudinal and transverse response functions from 3 He and 4 He is being 
conducted at the Mainz Microton (MAMI) facility by the Al-collaboration p|-|TT||. Their 
extraction of "experimental" spectral functions and of momentum distributions relies on 
a plane- wave- impulse-approximation (PWIA). In such an approximation the (e, e'p) cross 
section is proportional to the nucleon spectral function times an off-shell electron-proton 
cross section (cr ep ). Experimental analyses of this reaction employ, almost exclusively, the 
de Forest's ccl prescription for a ep with both nucleon form factors unmodified from their 
free-space form . 

Stimulated by this new experimental thrust, we report here relativistic plane- wave- 
impulse-approximation (RPWIA) calculations of the (e, e'p) cross section in the quasielastic 
region. Our motivation for such a study is fourfold. First, we employ an established RPWIA 
formalism, first introduced in Ref. [^] and recently extended to the kaon-photoproduction 
reaction |T4]Jl5l , for the study of the (e, e'p) reaction in the quasielastic region. Second, we 
use this formalism to compute the spectral function of 4 He in anticipation of the recently 
measured, but not yet fully analyzed, Al-collaboration data from Mainz |7HTll. Third, we 
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take advantage of the L/T separation at Mainz to introduce what we regard as the cleanest 
physical observable from which to extract the nucleon spectral function. Lastly, we high- 
light some of the challenges and opportunities that remain in the calculation of quasifree 
observables. 

There is a vast amount of literature on (e, e'p) reaction in the quasifree region. Most rel- 
evant to our present discussion is the one pertaining to fully relativistic calculations [[T~6| |28 



An extensive set of these relativistic studies has been conducted by the "Spanish" group 
of Udias and collaborators fl8|-f24"|| . These studies have shown that the many subtleties in- 
trinsic to the relativistic approach challenge much of the "conventional wisdom" developed 
within the nonrelativistic framework and that, as a result, a radical revision of ideas may 
be required. Relativistic effects originating from medium modifications to the lower compo- 
nents of the Dirac spinors and from the negative-energy part of the spectrum seem to play 
an important role in the quasifree process. Indeed, the much debated issue of short-range 
correlations at large missing momenta P§|-|3"Tf can now be attributed, at least in part, to 
contributions arising from the negative-energy states [|2l],[32 . 



The power of the theoretical approach employed here lies in its simplicity. Analytic 
expressions for the response of a mean-field ground state may be provided in the plane-wave 
limit. The added computational demands placed on such a formalism, relative to that from 
a free on-shell proton, are minimal. The formalism owes its simplicity to an algebraic trick, 



first introduced by Gardner and Piekarewicz [13], that enables one to define a "bound" (in 



direct analogy to the free) nucleon propagator. Indeed, the Dirac structure of the bound 
nucleon propagator is identical to that of the free Feynman propagator. As a consequence, 
the power of Feynman's trace techniques may be employed throughout the formalism. 

The paper has been organized as follows. In Sec. [Ij some of the central concepts and ideas 
of the semi-exclusive (e, e'p) reaction are reviewed. Special emphasis is placed on defining 
the bound-state propagator and the simplifications that this entails in the plane- wave limit. 
In Sec. [Til] we present our results for 4 He and discuss a (fairly) model-independent method 
for extracting the nucleon momentum distribution. Finally, a summary and conclusions are 
presented in Sec. [IV. 



II. FORMALISM 



In Refs. |14| , |15| a general formalism has been developed for the study of a variety of 
quasifree processes in the relativistic plane- wave impulse approximation (RPWIA). This 
formalism is now applied to the (e, e'p) reaction in a mean-field approximation to the Walecka 
model ||33|| . Although the use of a mean-field approach for a nucleus as small as 4 He is 
questionable, we allow ourselves this freedom in order to establish a baseline against which 
more sophisticated approaches may be compared. 

Following a standard procedure, an expression for the unpolarized differential cross sec- 
tion per target nucleon for the (e, e'p) reaction is derived. We obtain, 

/ d 5 a \ 4a 2 Ik 



P'll-Ml 2 - (1) 



dE' e dn v >dVt pl J lab Q 4 |k 



In the above expression k, k, and p' denote the linear momentum of the incoming electron, 
outgoing electron, and knocked-out proton, respectively. The four-momentum transfer is 
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defined in terms of the energy loss (u = E e —E' e ) and the three-momentum transfer (q = k— k') 
as Q 2 = q 2 — u 2 . The transition matrix element A4 is given in a relativistic mean-field picture 
by 



l-MI 2 



l» v W, 



(k'»k v + Wk' v - g" u (k ■ k' 



(2a) 
(2b) 



(2c) 



Here U(p', s') is the free Dirac spinor for the knocked-out proton, normalized according 
to the conventions of Bjorken and Drell while U am (p) is the Fourier transform of the 
relativistic spinor for the bound proton. Note that a denotes the collection of all quantum 
numbers necessary to specify the single-particle orbital, except for the magnetic quantum 
number (m) which is indicated explicitly. We have also introduced a "bound-state propaga- 
tor" 



w " v = 4 ( 2 / + i) ^(^ (p/ ' /)jAt ^ m(p) ) (^(pV)rz4m(p; 

= jT r ((f + M)rS a (p)f 



(3) 



normalized according to: 

d 3 p 



(27T) 



d 3 p 
(2vr) 3 



WL(P)*4 



1 . 



(4) 



Here j is the total angular momentum quantum number and 2j + 1 is the multiplicity of 
protons in the struck shell. It follows from simple kinematical arguments that the missing 
momentum p = p' — q is, in a mean-field picture, identical to the momentum of the struck 
proton. It is the possibility of mapping the nucleon momentum distribution that makes the 
(e, e'p) reaction so appealing. 

We now invoke an algebraic trick first introduced in Ref. JT3] to simplify the expression 
for the hadronic tensor W^ u . This technique is useful in quasifree processes as it enables 
one to cast the bound-state propagator of Eq. (^|) into a form identical in structure to that 
of the free Feynman propagator. That is, 

(5) 



S a (p) = {t a + M a 

where we have defined mass- and four-momentum-like \p% = (E a , p a )] quantities according 
to 



M a = - 



Pa = — 



glip) - fM 
gl(p) + fa(p) 

2ga(p)fa(p)p 



(6a) 

(6b) 

(6c) 
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Moreover, they satisfy the "on-shell relation" 



Pl = El-vl = Ml. (7) 
In these expressions g a (p) and f a (p) are the Fourier transforms of the upper and lower 



components of the bound-state Dirac spinor, respectively [IB ]. Using this form of the bound- 
state propagator the hadronic tensor simplifies to: 

= \T r ( (f + M) f (j a + M a )f). (8) 

The obvious similarity in structure between the free and bound propagators results in 
an enormous simplification: powerful trace techniques developed elsewhere may now be 
employed here to compute all (e, e'p) observables. Although the focus of this paper is the 
unpolarized cross section [Eq. (fj)] the formalism may be extended without difficulty to the 
case in which the electron, the outgoing proton, or both, are polarized. Yet, in order to au- 
tomate this straightforward but lengthy procedure, we rely on the FeynCalc 1.0 |35| package 



with Mathematica 2. to calculate all the necessary traces. For a general electromagnetic 
current operator for the proton, the output from these symbolic manipulations is trans- 
parent enough so that the sensitivity of the cross section to the various quantities in the 
problem may be assessed. Indeed, such a simplification will prove useful later in identify- 
ing the optimal observable from which to extract the spectral function. It is important to 
note, however, that this enormous simplification would have been lost had distortions been 
included in the formalism. Even so, the plane-wave approach discussed here, and used in 
most experimental extractions of the spectral function, is qualitatively useful. Moreover, if 
the main effect of distortions is to induce an overall suppression of the cross section without 
affecting significantly the distribution of strength, the plane-wave formalism provides solid 
quantitative predictions for a variety of spin observables fL4|,15 . 



Yet an important open question remains: what constitutes a suitable form for the nucleon 
electromagnetic current? A ubiquitous form given in the literature is 

/(?) = fi(?V + ^V F |. (9) 

While this form is certainly general, as only two form factors are required to fully specify the 
electromagnetic current for an on-shell nucleon, the form is not unique. Indeed, many other 
forms — all of them equivalent on-shell — may be used. For example, through a Gordon 
decomposition of the current one arrives at 

r(q) = (F 1 + F 2 )Y ~ Fi^T ■ (10) 

However, as soon as one of the nucleons goes off its mass shell, an off-shell choice must 
be made. This decision is crucial, as various on-shell equivalent choices may yield vastly 
different results. This off-shell ambiguity remains one of the most serious obstacles in the 
field. Several attempts have been made in the literature to overcome this hurdle. Perhaps 
the most celebrated treatment is due to de Forest who uses physical constraints, such as 
current conservation, to reduce this ambiguity |12|. He imposes this condition on the two 



forms of the electromagnetic current given above [Eqs. (§) and |10| and produces what are 
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known in the literature as the cc2 and the ccl forms, respectively. Although noteworthy, this 
effort does not resolve the ambiguity. For example, there is no unique way to impose current 
conservation; one may eliminate either the time component or the longitudinal component 
of the three- vector current [21]. Alternatively, one may adopt some guiding principle, such 
as vector-meson-dominance, to go off the mass shell. Here we adopt the "natural" choice by 
simply extrapolating off the mass shell the cc2-form, without imposing further constraints 
on the single- nucleon current. 



III. RESULTS 

As de Forest has done in the past, we now attempt to impose some approximate form of 
gauge invariance. Yet rather than concentrating on the nucleon current, we focus directly on 
the nuclear responses. First, however, we address some important issues in this regard. For 
any mean-field treatment of the (e, e'p) reaction to be gauge invariant, the mean-field poten- 
tial for the bound proton must be identical to the distorting potential for the emitted proton. 
This represents a challenging task. Indeed, mean-field approximations to the nuclear ground 
state give rise to real, local and energy-independent potentials that are in contradiction to 
the complex and energy-dependent potentials that are needed to describe the propagations 
of the outgoing proton. Thus, present-day calculations of (e, e'p) observables are presented 
with a dilemma. Calculations that use the same (real and energy-independent) mean fields 
to generate both the bound single-particle wave-function and the distorted wave satisfy 
gauge invariance but miss some of the important physics, such as absorption, that is known 
to be present in the outgoing channel. On the other hand, calculations that incorporate 
the correct physics via a phenomenological optical potential are known to violate current 
conservation . We offer here no solution to this complicated problem. Rather, we impose 



gauge invariance "ad-hoc" by adjusting the effective nucleon mass of the emitted proton so 
that the "gauge-variance" term, q^W^, be minimized. This procedure, with perhaps its 
unexpected outcome, is displayed in Fig. |I[ It shows that by decreasing the proton mass by 
about 20 MeV, one can restore gauge invariance in the calculation: q fl q u W tlL/ = 0. Although 
by no means fundamental, this "poor-man" distortion ensures the conservation of gauge 
invariance without compromising the clarity of the formalism. 

The essence of the experimental extraction of the spectral function is based on a nonrel- 
ativistic plane- wave result [36] : 



5 ^' p) P 'E' p a eN dE> e dn k ,dE> p dn p , ' 

However, this procedure is problematic. First, the quasifree cross section [the numerator 
in Eq. (|TT|)1 suffers from the off-shell ambiguity; different on-shell equivalent forms for the 
single-nucleon current yield different results. Second, the problem gets compounded by the 
use of an elementary electron-proton cross section (cr e jv) evaluated at off-shell kinematics jI2 



Finally, the projection of the bound-state wave-function into the negative-energy sector as 
well as other relativistic effects spoil the assumed factorization of the cross section derived 
in the nonrelativistic limit 21] . 
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Insights into the role of relativistic corrections, particularly those concerned with 
negative-energy states, may be gained by introducing the completeness relation in terms 
of free (plane-wave) spinors: 

£ [w(p, s) W(p, s) - V(p, s)V(p, s)]=l. (12) 

Naively, one would expect that the projection of a positive-energy bound state into a 
negative-energy plane-wave state would be vanishingly small. This, however, it is not the 
case ]32|] . At the very least one must recognize that the positive-energy plane- wave states, by 



themselves, are not complete. Moreover, it has been shown that the projection of the bound- 
state spinors into the negative-energy states dominate at large missing momenta and may 
mimic effects perceived as "exotic" from the nonrelativistic point of view, such as an asym- 



metry in the missing-momentum distribution [13] or short-range correlations [32]. Indeed, 



Caballero and collaborators have confirmed that these contributions can have significant 



effect on various observables, especially at large missing momenta [21 . 

To "resolve" the off-shell ambiguity it has become ubiquitous in the field to use the 
de Forest ccl prescription for evaluating the elementary cross section a e N — irrespective 
of the form of the electromagnetic current adopted to compute the quasifree cross section. 
This is the standard procedure used in comparing theoretical calculations of the spectral 
function to experiment. We may elect here to conform to tradition and use the de Forest 
ccl prescription to compute a e jq in Eq. (|TT|), but at a cost. A price must be paid because 
of the inconsistency in using one prescription for evaluating the single-nucleon current a e N 
and a different one (cc2) to evaluate the quasifree cross section. To illustrate this point we 
display in Fig. ^ the proton momentum distribution defined by 

p 2 (p) = Js(E,p)dE. (13) 

Note that the subscript "2" in p2 stands for two-body breakup. The graph displays the 
"canonical" momentum distribution (solid line) obtained from the Fourier transform of the 
IS 11 / 2 proton wave-function [see Eq. (|6"T5|)]. Note that this canonical momentum distribution 
has been normalized, as it is done experimentally, to the total number of protons in the shell 
(2 for the case 4 He). The other two curves were extracted from the quasifree cross section 
by adopting either the de Forest ccl choice for a e N (dashed line) or the cc2 prescription 
(dot-dashed line). In both cases the quasifree cross section has been computed using the 
"vector-tensor" form of the electromagnetic current, as given in Eq. (|H). The inset on the 
graph shows the integrand from which the occupancy of the shell may be computed. It 
is evident that the conventional ccl prescription of de Forest greatly overestimates p2 (it 
integrates to 3.6). We attribute this deficiency to the lack of consistency: the quasifree 
cross section has been evaluated using the cc2 form of the current while the elementary 
amplitude uses the ccl form. One can improve the situation by adopting the cc2 form in 
the evaluation of both. Yet significant differences remain; while the off-shell ambiguity has 
been reduced, it has not been fully eliminated. Moreover, the factorization assumption is 
only approximate, as it neglects the projection of the relativistic wave function onto the 
negative-energy spectrum and other relativistic effects. 

While a consistent relativistic treatment seems to have spoiled the factorization picture 
obtained from a nonrelativistic analysis, and with it the simple relation between the cross- 
section ratio and the spectral function [Eq. flOP I, the situation is not without remedy. Having 
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evaluated all matrix elements of the electromagnetic current analytically in the plane-wave 
limit, the source of the problem can be readily identified. Upon evaluating the coincidence 
cross section, one learns that the off-shell ambiguity is manifested in the form of several 
ambiguous "kinematical" factors. For example, one must decide what value to use for the 
energy of the struck proton. Should it be the binding-energy of the struck proton or should 
it be the on-shell value? This is not an easy question to answer. Energy conservation 
demands that the energy be equal to the binding energy (E^ — E'—u) yet the equivalence 
between the various forms of the electromagnetic current is derived assuming the on-shell 
dispersion relation (E p = y/p 2 + M 2 ). This is one of the many manifestations of the off-shell 
ambiguity: kinematical terms that are well defined for on-shell spinors become ambiguous 
off-shell. In Ref. [[[2] de Forest resolves the ambiguity, by fiat, using the on-shell choice. 



Perhaps a better option may be looking for an observable, that even though might be more 
difficult to isolate experimentally, it may display a weaker off-shell dependence than the 
unpolarized cross section. To do so we examine the various components of the hadronic 
tensor. We find, perhaps not surprisingly, that the longitudinal component of the hadronic 
tensor could be such a model-independent observable. Ignoring (for now) the anomalous 
part of the electromagnetic current, the Dirac-Dirac component of the longitudinal tensor 
[see Eq. (ffl)] becomes: 



oo 

DD 



F? \M a M -p a -p' + 2E a E' p ] = Ft \M a M + E a E' p + p a ■ p' 



(14) 



This expression depends exclusively on p a and p', which are unambiguous. Note that for 
scattering from a free on-shell nucleon the above expression becomes: 



tree 



M 2 - p ■ p' + 2E P E' =Ft M 2 + E P E'+ p • p' 



(15) 



Also note, as a consequence of the lower component of the bound-state spinor f a (p) being 
substantially smaller than the upper component g a (p), that |p Q | <^E a while M a ~E a . This 
is true even though the lower-to- upper ratio f a /g a has been enhanced considerably in the 
nuclear medium relative to its free-space value. This is an important step towards isolating 
an observable sensitive to the spectral function. Indeed, if the longitudinal component of 
the hadronic tensor is computed in parallel (p' = q) kinematics, Eqs. (14) and (|l5|) reduce 
to the following simple expressions: 



W 00 

^DD 



IV 00 



free 



F 2 (E' p + M) 



F 2 (E'+M) 



71 



2" p 



1 ± 



' UpY 



IP 



(E p +M) 



1 ± 



iPl 



IP 



E p + MJ \E p + M / 



(16a) 
(16b) 



The ± sign in the above expressions corresponds to a missing momentum p either parallel 
or antiparallel to p'. We observe that up to second-order corrections in the small (lower-to- 
upper) ratios, the hadronic tensor is proportional to the energy-like (or mass-like) quantity 
given m Eqs (|BD. Yet this energy-like quantity E a is nothing but the Fourier transform 
of the bound-state nucleon density. Thus we conclude that, in a mean- field treatment, the 
nucleon spectral function is proportional to the longitudinal response. That is, S(E,p) oc 
W^ocEa. Thus, the (Dirac-Dirac component of the) longitudinal hadronic tensor is, up to 



8 



second-order corrections in the lower-to-upper ratios, proportional to the nucleon spectral 
function. Indeed, the nucleon momentum distribution may now be easily extracted from the 
longitudinal response. It becomes 



p 2 = 2(2j + l) (E P + M)(W°° D /W, 



oo 

DD 



free 



(17) 



The momentum distribution for 4 He is displayed in Fig. ^| using various methods for its 
extraction. The solid line gives the "canonical" momentum distribution, obtained from 
the Fourier transform of the IS* 1 ' 2 proton wave-function [see Eq. fl6T5|)1. The momentum 
distribution extracted from the longitudinal response as defined in Eq. (|T7|) (dot-dashed line) 
is practically indistinguishable from the canonical momentum distribution. While it appears 
that a suitable observable has been found from which to extract the nucleon momentum 
distribution, it may be argued, and justifiably so, that Wjj D is not a physical observable (as 
F 2 has been neglected). Hence, the merit of such an extraction may be put into question. 
To show that the above procedure is still robust, we display in the figure (with a dotted line) 
the momentum distribution extracted from the full longitudinal response, namely, one that 
also includes the anomalous component of the current. This result remains indistinguishable 
from the canonical momentum distribution. Although this behavior is general, it is most 
easily understood by limiting the discussion to the case of parallel kinematics. In this case 



the longitudinal response becomes equal to [13 



R h = W m = (E' p + M) 



71 



(F 1 -i' p qF 2 )±(i' p F l + qF 2 



' Up) ' 



(18) 



The contribution from the anomalous form factor F 2 to the longitudinal response is small 
because it appears multiplied by two out of three "small" quantities in the problem: the lower- 
to- upper ratio, £p= \ p'\/(E' + M), and q= |q|/2M. Thus, up to second order corrections in 
these small quantities, the longitudinal response is given by 



R L ~F?(E' p + M) 



71 



P 



Ft(E' p + M)E a 



(19) 



The last calculation displayed in Fig. [| corresponds to a momentum distribution extracted 
from the factorization approximation using the cc2 form for the electromagnetic current 
(long dashed line). The momentum distribution extracted in this manner overestimates the 
canonical momentum distribution over the whole range of missing momenta and integrates 
to 2.9 rather than 2; this represents a discrepancy of 45 percent. 

In summary, the longitudinal response appears to be a robust observable from which to 
extract the nucleon momentum distribution. Experimentally, one should proceed as follows: 
perform a Rosenbluth separation of the (e, e'p) cross section so that the longitudinal response 
(Rl = W 00 ) may be extracted. This expression should then be divided by the corresponding 
single-nucleon response. Up to a simple and unambiguous kinematical factor this yield, at 
least in the plane- wave limit, the nucleon momentum distribution: 



p 2 = 2 (2j + 1) (E p + M) 



R 



free 



(20) 
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Note that up to second order corrections in various small quantities, this form is indepen- 
dent of the small components of the Dirac spinors and also of the negative-energy states. 
Moreover, it is also free of off-shell ambiguities. Indeed, we could have used the ccl form 
of the electromagnetic current and the results would have remained unchanged. We regard 
the outlined procedure as much more robust than the conventional one given in Eq. ( |TT|) 
because the transverse component of the hadronic tensor is strongly dependent on the small 
components of the wave-function and also sensitive to off-shell extrapolations 



In Fig. £| a comparison is made between our results and nonrelativistic state-of-the-art 
calculations of the momentum distribution of 4 He. The solid line displays, exactly as in 
Fig. [I the canonical momentum distribution. We see no need to include the momentum 
distribution extracted from the longitudinal response [Eq. (|20|) 1 as it has been shown to 
give identical results. In addition to our own calculation, we have also included the varia- 
tional results of Schiavilla and collaborators |37j , for both the Urbana (dashed line) and 
the Argonne (long-dashed line) potentials, with both of them using Model VII for the 
three-nucleon interaction. The variational calculation of Wiringa and collaborators [f40| -|42|| 
(dashed-dotted) has also been included; this uses the Argonne vl8 potential |G| supple- 
mented with the Urbana IX three-nucleon interaction ||44|| . Figure |] also shows NIKHEF 
data by van den Brand and collaborators as well as preliminary data from MAINZ 

by Florizone and collaborators |F]|| for three different kinematical settings. (Results in final 
form will be submitted shortly.) Comparisons to the preliminary Mainz data of Kozlov and 
collaborators IPHTTfl have also been made (although the data is not shown). These measure- 
ments are consistent, in the region where comparisons are possible, to the experimental data 
of both van den Brand and Florizone. Thus, high-quality data for the momentum distri- 
bution of 4 He is now available up to a missing momentum of about 200 MeV. We find the 
results of Fig. ^ quite remarkable. It appears that a simple relativistic mean-field calculation 
of the momentum distribution rivals — and in some cases surpasses — some of the most 
sophisticated nonrelativistic predictions. The mean-field calculations reported here, with 
the scalar mass adjusted to reproduce the root-mean-square charge radius of 4 He, provide a 
good description of the experimental data. Still, theoretical predictions of the momentum 
distribution overestimate the experimental data by up to 50-60%. Part of the discrepancy 
is attributed to distortion effects which are estimated at about 12% lf7|,|47|. However, dis- 
tortions are not able to account for the full discrepancy. We have argued earlier that an 
additional source of error may arise from the factorization approximation [see Eq. (|ll])] used 
to extract the spectral function from the experimental cross section. The use of an off-shell 
prescription, such as the ccl prescription for a e N, combined with the in-medium changes in 
the lower-component of the Dirac spinors contaminate the extraction of the spectral func- 
tion. One could estimate the source of the off-shell ambiguity by monitoring the variations 
in the spectral function as other on-shell equivalent forms for the single-nucleon current 
are used. While such an approach is useful for estimating a theoretical error, it is clearly 
not sufficient to eliminate it. We are confident that the approach suggested here, based 
on the extraction of the spectral function from the longitudinal response, is robust. While 
the method adds further experimental demands, as a Rosenbluth separation of the cross 
section is now required, the extracted spectral function appears to be weakly dependent on 
off-shell extrapolations and relativistic effects. If deviations between experiment and theory 
still persist, these may suggest physics beyond the baseline model, such as violations to the 
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impulse approximation or to the independent particle picture. 



IV. CONCLUSIONS 

To summarize, we have calculated the spectral function of 4 He in a plane-wave approxi- 
mation to the (e, e'p) reaction using a fully relativistic formalism. We have taken advantage 
of an algebraic trick originally introduced by Gardner and Piekarewicz and of our recently 
developed relativistic formalism for quasifree processes to arrive at transparent, analytical 
results for the quasifree reaction. We have found that a simple relativistic mean-field calcu- 
lation of the momentum distribution in 4 He rivals — and in some cases surpasses — some 
of the most sophisticated nonrelativistic predictions to date. These calculations attempt to 
provide theoretical support to the recently measured, but not yet fully analyzed, Al col- 
laboration data from Mainz. The final experimental reports are expected to be published 
shortly. 

We have also demonstrated that a more robust procedure, relative to the conventional 
factorization prescription, exists for extracting the spectral function. This procedure uses 
the ratio of quasifree to single-nucleon longitudinal responses, rather than the ratio of cross 
sections, to isolate the momentum distribution. We have shown that the longitudinal ratio 
is fairly insensitive to off-shell ambiguities and to the negative-energy part of the spectrum, 
as both of these effects appear as second-order corrections to a "canonical" momentum 
distribution. This ceases to be true in the case of the ratio of cross sections because the 
transverse response is sensitive to both effects. While this procedure relies on a Rosenbluth 
(L/T) separation of the quasifree cross section, and thus presents the experimentalist with 
a more demanding task, the experimental field has evolved to such a level of maturity that 
L/T separations are now almost routine. Indeed, in a recent publication |48|] a Rosenbluth 



separation of the 3 He(e, e'p) cross sections was made in order to extract "longitudinal" and 
"transverse" spectral functions in the hope of resolving the anomaly in the longitudinal- 
transverse ratio alluded to in the introduction. We speculate that the sensitivity of the 
transverse response to more complicated dynamical processes might be partially responsible 
for the quenching of the longitudinal-transverse ratio. 

Finally, although in this article we focused exclusively on the spectral function, the 
formalism presented here may be extended in a straightforward fashion to the calculation of 
spin observables in quasifree electroproduction processes. Indeed, we speculate that, because 
the ratio of quasifree cross sections are fairly insensitive to distortion effects, spin-observables 
may be a more fruitful testing ground for our relativistic plane-wave model. Moreover, our 
formalism may be easily extended to neutrino-induced reactions. It has been suggested that 
a measurement of the ratio of neutral to charge- changing neutrino-nucleon scattering may 



provide a clean signature of the strange-quark content of the nucleon ||9| . This measurement 
is believed to be free from most of the uncertainties, such as radiative corrections, that hinder 
the parity-violating electron scattering program. Yet neutrino experiments suffer from very 
low counting rates. To remedy this situation neutrino experiments employ large quantities of 
nuclear targets (such as organic scintillators) that provide both the target and the detection 
medium. Thus neutrinos interact, not only with the free protons in the target, but also 
with protons and neutrons bound to nuclei; hence, one must compute quasifree (u, v' p) and 
(v, yTp) cross sections. (Of course, one must integrate the quasifree cross section over the 
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undetected outgoing neutrino). Therefore, the relativistic plane- wave formalism presented 
here is ideally suited, after including an additional axial-vector term in the single-nucleon 
current, to predict ratios of quasifree neutrino- nucleus cross sections in the quasifree region. 
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FIG. 1. The gauge variance term q ll q u W tlu for 4 He as a function of the scattered proton mass 
calculated in parallel kinematics for an incident photon energy of E- mc = 855 MeV and a momentum 
transfer of q = 685 MeV. The right panel is a magnification of the boxed area in the left panel. 
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FIG. 2. The proton momentum distribution p2 for 4 He as a function of the missing momen- 
tum calculated at an incident photon energy of E[ nc = 855 MeV and a momentum transfer of 
q = 685 MeV. The solid line is the relativistic mean-field calculation, while the dashed and 
dot-dashed lines display the momentum distribution extracted from a factorization approxima- 
tion using the ccl and cc2 prescriptions for cr e N, respectively. The inset shows the corresponding 
integrands from which the shell occupancy may be extracted. 
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FIG. 3. The proton momentum distribution p2 for 4 He as a function of the missing momen- 
tum calculated at an incident photon energy of E[ nc = 855 MeV and a momentum transfer of 
q = 685 MeV. The solid line is the relativistic mean-field calculation, while the dashed and 
dot-dashed lines display the momentum distribution extracted from the longitudinal response Rl 
without including (dot-dashed) and including (dotted) the contribution from the anomalous form 
factor F2. Finally, the dashed curve is obtained by using the factorization approximation with 
the cc2 prescription for a e M- The inset shows the corresponding integrands from which the shell 
occupancy may be extracted. 
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FIG. 4. A comparison between our relativistic calculations, nonrelativistic calculations reported 
elsewhere, and experimental data for the proton momentum distribution in 4 He. The solid line is 
our mean-field calculation while the dotted curve is our calculation using the factorization approx- 
imation at incident photon energy of E- inc = 855 MeV and a momentum transfer of q = 685 MeV. 
The nonrelativistic calculations of Schiavilla et al, are included for both the Urbana (dashed) and 
the Argonne (long dashed) potentials as well as the calculations of Wiringa et al, (dashed-dotted). 
The NIKHEF data of van den Brand et al, for two different kinematical settings as well as pre- 
liminary data of Florizone et al, (Al collaboration) which were measured at MAMI (Mainz) are 
also shown. 
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